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I discuss some of the forces acting on vortices in charged superfiuids, paying particular 
attention to the way that the Berry and Aharonov-Casher phases combine to reflect the 
classical magnetohydrodynamics. 

1 Introduction 

The motion of vortices is the only significant source of dissipation in a supercon- 
ductor, making it a topic of some technological importance. Much of the modern 
work on vortex dynamics relates to many-vortex interaction effects so it is rather 
surprising to find that there is still considerable debate about the forces acting on 
an individuaLvprtex. This debate was sparked by the claim in papers by Ao, Thou- 
less and Niulila that there exists a universal, exact expression for the transverse 
Magnus force on a vortex. Arising from topological effects, the magnitude of this 
Magnus force should not be affected by impurities or quasiparticles except in that 
they change the value of the superfluid density. 

ThejA^IhouJess-Niu papers stimulated a number of authors to re-examine the 
problemaui3QQi3'Ll with conclusions ranging between support for their viewpoint, 
through the suggestion that spectral flow, yet another topological effect, partly 
cancels the Magnus force, to claims that there are no topological effects at all. 
It has-.also led to an ingeneous experiment aimed at a direct measurement of the 
force£9. 

Given the interests of the audience at this conference, I will focus on only one 
aspect of this problem — the effect, if any, electromagnetic interactions have on 
the topological Magnus force. In the next section I will review the connection 
between the topological phase and the Magnus force for neutral superfiuids, and 
then show how the phase accounting is modified by Aharonov-Casher phases due to 
the motion of the magnetic flux which is tied to the vortices. In the third section I 
will show how the various phase cancellations reflect a purely classical re-routing of 
the momentum flux from the charged superfluid to the background lattice. A final 
discussion section will, for completeness, briefly review how spectral flow induces 
relaxation effects that need to be incuded to complete the picture. 

2 Berry and Aharonov-Casher phases 
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2.1 Neutral Fluids 



Variation of the action 



S = J dtd 3 x j-i(^¥> - + 2^|V^| 2 + ^(M 2 - po) 2 } (1) 

gives rise to a non-linear Schrodinger equation 

td t <p = - 7 }-V 2 <p + \(\p\ 2 - Po )cp (2) 
2m 

describing fluctuations of the scalar field ip about a stationary configuration with 
\(p\ 2 = pa. This equation is called the Gross-Pitaevskii (GP) equationO and is 
often usecLas a simple model for the motion of a Bose condensate. Long ago 
Madelunglj showed that the real and imaginary parts of the Schrodinger equation 
can be interpreted as the equations of mass and momentum conservation for a simple 
fluid, so it is no surprise that vortex solutions to the GP equation are consistent 
with Kelvin's circulation theorem from elementary fluid dynamics — i.e the vortices 
are advected along with the flow. For example, in three dimensions, there exists 
a family of smoke-ring-like vortex solutions to- the GP equation which move along 
their symmetry axis without change of radiust2l. Any attempt to force a vortex to 
move with respect to the flow requires the application of a transverse lift force— .In 
the superconductivity literature this is usually referred to as the Magnus forcet-3. 

We can simultaneously obscure and illuminate the existence of this force by 
some cosmetic rewriting of the action. We begin setting cp — ~Jpe lB to get 



S = j dtd 3 x^pd t 9 + £-(X79) 2 + ^(p-p ) 2 + ...^ 



(3) 



(where the dots indicate a presently uninteresting term depending on the gradient 
of p.). Next we integrate over p to find (again ignoring the gradient of p terms) 



S = J dtd 3 x ^ Po d t + 7^(y0) 2 - ^(d t 9 + \mv 2 s ) 2 ^ (4) 

where v s = ^V6* is the supefluid velocity. 

The first term is a total derivative and therefore does not contribute to the 
classical equation of motion for 9. If we temporarily suppress the v 2 part of the 
third term, variation of the action gives rise to a wave equation which describes 
sound waves propagating in the fluid. We should not, however, use this truncated 
equation for the vortex dynamics. Wave-equation dynamics for an angular variable 
like 9 appears in the relativistic form of the X — Y model and it. is well known that in 
such a model a vortex loop will rapidly shrink and disappeared. Throwing awajz-jthe 
higher order terms in v s , although common in the superconductivity literatural3, is 
therefore not safe when one is interestetLin vortex motion. This because omitting 
these terms destroys galilean invarianceO. 

It is not safe to discard the first term pod t 9 eitherTf one wishes to use S to 
compute the partition function via an imaginary time path integral over configura- 
tions with periodicity (3 = 1/kT, the first term will give rise to a topological phase. 

This is done, for example, in the otherwise excellent texbook by Popov Q. 
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Although the field <p must return to its original value at r = f3, the angle 9 will in 
general have wound through 2ttN where N is an integer. In two dimensions the or- 
der parameter phase associated with a vortex contains a part 6(r, t) = Arg(r— vi,{t)) 
therefore moving the vortex in a closed path will contribute a phase factor equal to 

exp{2iripo(Area)} (5) 

where {Area) is the area enclosed by the vortex trajectory. This is exactly the 
same phase factor that would be generated by a charged particle with interaction 
action <f> ex^A^ dt moving in a uniform magnetic field. It must therefore signal 
the existence a transverse "Lorentz" force of magnitude 2irpoV vor t ex . Taking into 
account that the circulation around the vortex is k — 2n/m this is exactly the 
expected magnitude of the classical Magnus force p§mKV vort&x . 

This connection between this topological "Berry" phase and the Magnus force 
in a superfluid was pointed out in Rcf.l. It became the motivation for a more 
general theorem on the existence of transverse forces on vorticesB. 



2.2 Charged Fluids 

What happens to the vortex-induced Berry phase when the fluid is composed of 
charged particles? To find out it is simplest to perform a duality map which raises 
the status of the vortices from mere defects in the order parameter to the central 
objects of discussion. There are many references for such maps, but since the 
notation and inspiration for this section comes directly from a papenla written by 
a member of the audience I will recommend that for an introduction. 

For geometric simplicity I will restrict myself to two dimensions. In the following 
roman subscripts such as a will run over the two space dimensions while greek 
subscripts such as a will run over both euclidean time (a — 0) and space dimensions. 

When we couple our scalar field to a gauge field the Lagrange density becomes 

L = i(^(ao-i e Ao)0-0(9o + ^ o )^) + ^|(9 a -i e AJ0| 2 + ^(|^| 2 -po) 2 + ^ 2 . 

(6) 

Set (f = fe t9 to write this as 

L = if 2 (d 6 - eA a ) + ^(d a f 2 ) 2 + ^(d a 9 - eA a f + ±{f - Pa ) 2 + -F 2 . (7) 
Now introduce a pair of Hubbard-Stratovich fields C a and also set Co = f 2 to get 

k A v - eAJ + J- Wg£ + + \ ( c. - ,f + (8) 

Next we isolate the vortex part of the phase 9 by writing 9 = 6 + r\ where 
9 = J2i Arg(r — r^)) is the singular part of the phase due to vortices at r,, while rj 
is the remaining non-singular part. Integration over T] enforces the mass conservation 
equation 8 t Cq + d a C a = 0. This is automatically satisfied if we write 

C M = epvadubo-. (9) 
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Regarding b a as a gauge potential then suggests renaming 

C a E a Co = H, (10) 
a set of pseudo electric and magnetic fields. After defining 

= e^dudaO (11) 
to be the vortex 3-current density we find that the Lagrange density takes the form 

L' = ib„(K„ - ee^dvAa) + ^\E\ 2 + -{H - Po ) 2 + ^F 2 + ... (12) 

where the dots represent the same higher order gradients of the density that were 
omitted in the earlier section. Assuming p ~ po, we can write the Lagrange density 




where F* = \e^ va F va is the dual of the real electromagnetic flux. 

We see that the vortex current is coupled to a gauge field E that apparently 
mediates coulomb interactions between vortices but actually accounts for the at- 
tractive and repulsive forces due to the kinetic energy of the fluid, and to a magnetic 
field H which has an equilibrium value po even when the vortex is at rest. Motion 
with respect to this background field gives rise to a Lorentz force which is just 
the Magnus force in this pseudo-field language. We also see that the source for 
the fields is not the bare vortex current but the vortex current less the flux of the 
genuine magnetic field that acompanies the vortex core. As in any Abrikosov flux 
tube solution, this magnetic flux completely screens the vortex current and appears 
to eliminate both the "coulomb" interaction between the vortcies and the Magnus 
force. The coulomb screening is a geniunc effect - the dual of the Anderson-Higgs 
mechanism. Is the Magnus force banished as well? The answer is no! We have 
omitted an important term in Eq.[j| This is the interaction of the electromagetic 
field with the background ions. Consideration of this is necessary if we are not to 
have infinite coulomb energy. Including an Li ons = iepi on Ao of the standard form 
and writing the ionic 3-current (a static charge of course) as the curl of a field b l ° n 
the final form of the vortex field interaction becomes 




When the fluid is at rest with respect to the background charges away from the 
vortices (which it has to be unless one is within a penetration depth of the boundary 
of the system) the two terms prortional to F* cancel and the Magnus force is 
restored. 

In terms of topological phases, what is happening is that the Berry phase from 
the superfluid is being supplemented by two Aharonov-Casher phasesEJ arising from 
moving the magnetic flux line round the charges in the system. One Aharonov- 
Casher phase comes from moving the flux through the mobile fluid charges and one 
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from moving it through the static background ion charges. Coulomb interactions 
force the net charge inside any macroscopic region to be zero so the two Aharonov- 
Casher phases must cancel. n 

The results of this section differ from those in some recent papersa. Although 
in the end they obtain the same total force on the vortex, the authors of these 
papers claim that inclusion of electromagnetic interactions causes the topological 
phase terms in the effective action to cancel. This cancellation originates in the 
form they chose for the interaction with the background charge. I have chosen to 
represent it as iepi on A^ whereas these papers argue that gauge invariance requires 
us to write this interaction term as a gauge covariant derivative iepi on (Ao — do9), 
where 9 is the phase of the electron order parameter. Then, because pi on = pa 
the terms in 8q9 cancel. I would argue that since the background ion current 
is conserved separately from the electron current, the conventional form is gauge 
invariant without ami. additional term. Similar opinions have been expressed by Zhu, 
Gaitan and VolovikEil. In the next section I will argue that my choice of interaction 
is the correct one by relating the phase cancellations to the actual mechanical forces 
acting on the system. 

3 Kutta-Joukowski Theorem for Charged Fluids 

The original example@ exhibiting the Aharonov-Chasher phase was constructed to 
be similar to the Bohm-Aharonov effect in that there was a quantum effect even 
though the classical force acting on the moving object vanished. However, just as 
the Bohm-Aharonov phase in a region with non vanishing magnetic field becomes 
path dependent and signals the existence of the Lorentz force, so the Aharonov- 
Casher phase in a region of non-vanishing charge density becomes path dependent 
and signals that moving flux through a charge distribution requires application of 
transverse force — the reaction force to the Lorentz force the charges feel in the 
moving magnetic field. (More accurately the charges feel the electric field arising 
from the Lorentz-transformed magnetic field). We can therefore avoid all contro- 
versy over cancelling phases and gauge invariance by focussing on the equivalent 
classical forces experienced by the fluid and the background ions. 

3. 1 Neutral Fluids 

To understand the forces on the vortex, or equivalently the momentum flux into and 
out of the vortex core, we may as well look at the forces on a hypothetical solid body 
occupyingj-the same region as the core. We will begin by reviewing the traditional 
derivations of the lift force on an aerofoil with circulation k which is being held 
stationary in electrically neutral, two-dimensional, incompressible potential flow. 
This force is most easily obtained from the momentum flux tensor. 

Tij = pmViVj + g^P. (15) 

Since we are interested in steady, irrotational motion with constant density we 
may use Bernoulli's theorem, P+^pm\v\ 2 — const., to substitute — ^pm\v\ 2 in place 
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of P. (the constant will not affect the momentum flow). With this substitution T^- 
becomes a traceless symmetric tensor 

Tij = pm{v t v 3 - -gij\v\ 2 ). (16) 

In two dimensions it is convenient to use complex coordinates z = x + iy where 
the metric tensor becomes g zz = g-zz — 0, g zz = g- z = \ and g zz — g zz — 2 etc. In 
these coordinates v z = ^(v x — iv y ) and 

T = T ZZ = -(T wx ~ T yy - 2iT xy ) = pm(v z ) 2 . (17) 

T is the only component of we will need to consider. T— is simply T while 
T z - — = T zz because Ty is traceless. 

Any body-force fi acting on the fluid is given by the divergence of the momen- 
tum flux tensor 

fi = 9 kl d k Tu (18) 

or, in our coordinates, 

U = 9"dgT„ + g z ~d z T zz = 2^T. (19) 

From this we see that, in a steady flow, the momentum flux P out of a region 
must be given by 

P z = — [ f z dzdz = - [ d^Tdzdz = - I Tdz. (20) 
2« J i J i Jan 

When there is no body-force, T is analytic and the integral will be independent of 
the choice of contour dil. 

To apply this result to our aerofoil we should take dfl to be its boundary. Then 
P z becomes minus the integral of the pressure over the boundary of the body, i.e the 
force exerted on the fluid by the aerofoil. Evaluating the integral is not immediately 
possible because the velocity V z will be a complicated function of the shape of 
the body. We can, however, exploit the contour independence of the integral and 
evaluate the integral over a contour encircling the aerofoil at large distance where 
the flow field takes the form 

v, = U, + ~ + 0(\). (21) 
4m z z^ 

To confirm that this flow has the correct circulation we compute 

v • dr = <j> v z dz + <j> ifjdz = k. (22) 

Subsituting in Eq.|20| we find that 

P z = —ipmnUz. (23) 



6 



In conventional coordinates the reaction force on the body F = P is 

F x = pmnUy 

F y = - pm,KU x . (24) 

The body therefore exerts a transverse lift force on the fluid proportional to the 
product of the circulation with the asymptotic velocity. This is the Kutta- Joukowski 
theorem. 



3.2 Charged Fluids 

How is this result modified when the neutral fluid is replaced by one composed 
of charged particles? We wish to model a charged superfluid where the superfluid 
velocity is obtained from the order parameter phase, 0, via the equation 

v=— (V0-eA). (25) 
to 

The flow is therefore no longer irrotational but instead has vorticity 

lo = V A v = -— V A A = -— B. (26) 
m m 

The relation eB + muj = gives rise to the Meissner effect: using the Maxwell 
equation V A H = J with J = pev we find 

V 2 B + /j — B = 0. (27) 

TO 

Eq.^ implies that both the magnetic flux and the vorticity decay exponentially 
away from a flux-creating disturbance. A vortex defect is such a disturbance. The 
2ir phase winding produces a solenoid-like current near the vortex core. This current 
creates a magnetic field which in turn induces a local vorticity that tends to screen 
that of the defect over the penetration length Ao = y 'm / '(pe 2 'po) . The circulation 
about the vortex is therefore path-dependent and tends to zero as the integration 
path becomes further from the vortex core. 

From the fact that the induced vorticity exactly screens the circulation round 
the bare vortex we can find the net flux 



/ TO / , , TO 27T 27T 

$ = / B 3 dxdij = / ujdxdy = = — . (28) 

el e to e 



For the Cooper pairs of a BCS superfluid m = 2x (mass of electron) and e = 2x 
(charge of electron), so this flux quantum is half that appearing in the quantum 
Hall effect. 

As before, we must not forget the charge of the background ions. Without 
a background charge pi on w p there would be a huge Coulomb field. Indeed any 
deviation from equality between the two charges will tend to zero over a Debye 
screening length so we will assume that pi on = p everywhere. Further any tendency 
of the superfluid to flow relative to this background charge will create a magnetic 
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field and the resulting Meissner effect will restrict the flow to within a penetration 
length Ao of the boundary of the system. To avoid this current inhomogeneity and 
to obtain a theorem analagous to that of Kutta- Joukowski we will therefore assume 
that the asymptotic fluid flow U is equal to that of the rigid motion background 
ions. This means that we are working in the rest frame of the vortex which is itself 
moving at a steady velocity — U with respect to both the asymptotic fluid and the 
background ion crystal. We also assume that our body f2 contains static charge 
with density equal to that of the surrounding fluid. 

Despite the non-zero vorticity it is easy to see that the momentum flux tensor 
remains = pm(viVj — hg^ \v\ 2 ). The momentum supplied by the aerofoil is 
therefore still given by 

P z = -<f Tdz, (29) 
i Jan 

but now T is no longer analytic. Indeed we have that 

dgT=±epBv s . (30) 

We must therefore take dfl as the boundary of the aerofoil. (From now on we write 
B for S3 since this is the only non-zero component of the field). Although we cannot 
use analyticity to send the integration contour off to infinity, we can still integrate 
by parts and find 

~ / Tdz = --[ d^Tdzdz (31) 
1 Jan 1 Jn 2 \n 



ep 
2 



Bv z dzdz. (32) 

R 2 \f2 



This last expression is simply the Lorentz force on the fluid outside the aerofoil. 
We now note that the Maxwell equation W/\H = J is in our complex coordinates 

id z B = ep p(v z — U z ) in R 2 \S! 

= —epopU z in f2. (33) 



So 



P z = — / {Bd z B)dzdz~— BU z dzdz. (34) 



The last term in Eqj34|is independent of the flow field outside the body but still 
depends on details of the system through the domain of integration. We will see 
that we can combine both terms to get an expression for the lift that is completely 
independent of any such details. 

We observe that 



— / (Bd z B)dzdz = — / d z -B 2 dzdz 

B 2 dz 



2 Mo Jr2\o 2^o Jn 2 \n 

2/^0 Jan 2 
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/ dJ-B 2 dzdz 



Jn 



— I epBU z dzdz. 
2 Jn 



(35) 



Putting the two parts together then gives 




BU Z dzdz = —iep$U z 



(36) 



This expression for the lift force on the aerofoil depends only on the total flux and 
the asymptotic current. Taking into account the relationship between the neutral- 
fluid quantum of circulation and the flux in the Abrikosov vortex we see that the 
force on the aerofoil, or equivalently the momentum flux into or out of a vortex, is 
unchanged by the introduction of electromagnetic interactions. 

What is happening is that in the neutral case momentum is flowing out of the 
vortex and increasing the net momentum of the fluid. In the charged case the 
moving vortex supplies exactly the same amount of momentum to the fluid, but 
the magnetic field acts as a clutch coupling the fluid to the background ions. Any 
relative motion creates a magnetic field that tends to inhibit it. The momentum 
is therefore transfered to the background ions. That is why the force the aerofoil 
exerts on the fluid is exactly equal to the Lorentz force the magnetic field exerts 
on the background ions as a result of their being dragged through the vortex's flux 
tube. A similar mechanism occurs in the quantum Hall effect where motion of 
either a Laughlin quasiparticle or a skyrmion-transfers momentum to the electrons, 
and they in turn hand it on to the magnetcJ. An alternative way of thinking of 
the same physics is to follow Nozieres and VineriO and realise that in the charged 
fluid the circulation, and hence the Magnus force, is reduced as we look at contours 
at larger and larger distance from the vortex. The lost Magnus force is however 
exactly compensated by the Lorentz force on the fluid within the contours. This 
re-routing of the momentum exactly mirrors the phase cancellations described in 
section 2. 

4 Discussion 

We have treated the superfluid as if it were a Bose condensate of Cooper pairs and 
therefore modeled its motion by the GP equation. This model differs in several ways 
from BCS electrons in a real metal. Firstly the GP equation is galilean invariant. 
This, in my opinion, is not a problem. One often models electrons in a metal by 
taking a parabolic band dispersion curve e = k 2 /2m*, and the resulting pseudo- 
galilean invariance provides a useful check on one's computation. 

Secondly some people object to modeling the dynamics by a first-order equation 
at all. For the_neutial case this should not be an issue. It has been shown by 
various authorsE-ffij'EEl that this is the correct form of the low energy dynamics for 
a galilean invariant BCS system at T = 0. The quantity tp appearing in the GP 
equation is not, however, the BCS order parameter. It is instead the combination 
<p = ^fpe %e where 9 is the phase of the order parameter and p is the fluid density. 
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(You may have noticed that throughout this talk I have never refered to the scalar 
field <p as the "order parameter" . I have only refered its phase as the "phase of 
the order parameter"!) These variables are those associated with the low-energy 
Goldstone mode. The magnitude of the BCS order parameter |A| is not associated 
with Goldstone-mode dynamics and so does not appear in the low-energy effective 
action. 

Finally, and most importantly BCS vortices have internal dynamics associated 
with quasiparticle core states. Volovik suggested that there should be a novel axial- 
anomaly driven spectral flow among these states in a moving vortex, and that 
this spectral flow would cancel most of the Magnus forcetl. It turns out that once 
one takes the discreteness of the core states into, account then one must include 
relaxation processes in order to see their effecl£3ll With relaxation included this 
turns_put to be a previously known effect known as the Kopnin-Kravtsov (KK) 
forceo, a sort of transverse friction which exhanges momentum with the crystal 
lattice by scattering of core quasiparticles off lattice defects. How does this extra 
transverse force fit in with the claims of Ao et al. that the Magnus force is all there 
is? The theorem of Ao-Thouless and Niua states that if a vortex moves against the 
flow, then momentum enters the superfluid at exactly the rate given by the Magnus 
force. When we include the KK friction a vortex placed in a flow field will find itself 
moving at a velocity different from the asymptotic flow. -Momentum must therefore 
enter the fluid at exactly the rate given by the theoremla. In order to solve for the 
relative speed of the vortex and the fluid, though, we need another expression for 
the rate at which momentum is entering the system. The actual mechanism by 
which the momentum is being gained is that momemtum is being transfered to the 
core states from the lattice by KK scattering, and from there flows out into the fluid 
in accord with Eq. |2^. The rate of KK scattering depends on the relative velocity 
of the vortex and the crystal lattice, so equating these two expressions for the rate 
of momentum transfer determines the steady state vortex velocity in terms of the 
relative velocity of the lattice and the flow field. 
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